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ABSTRACT 

Aims. The statistical properties of radio halos can be used to discriminate among the possible models for their origin. Therefore an 
unbiased and exhaustive investigation in this direction is crucial. 

Methods. With this goal in mind in this paper we revise the occurrence of radio halos in the redshift range 0-0.4, combining the low 
redshift (z < 0.2) statistical study of XBACs clusters with the NVSS (by Giovannini et al. 1999) with our recent results from the radio 
follow up of REFLEX and eBCS clusters, the GMRT radio halo survey, at higher redshift (0.2 < z < 0.4). 

Results. We find a significant statistical evidence (at 3.7cr) of an increase of the fraction of clusters with Radio Halos with the X-ray 
luminosity (mass) of the parent clusters and show that this increase is in line with statistical calculations based on the re-acceleration 
scenario. We argue that a fundamental expectation of this scenario is that the probability to have radio halos emitting at hundred MHz 
is larger than that at GHz frequencies and thus that future radio interferometers operating at low frequencies, such as LOFAR and 
LWA, should detect a larger number of radio halos with respect to that caught by present GHz observations. We also show that the 
expected increase of the fraction of clusters with radio halos with the cluster mass as measured with future LOFAR and LWA surveys 
should be less strong than that in present surveys. 
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. 1. Introduction 

■ Radio halos (RH elsewhere) are diffuse Mpc-scale radio sources 

■ observed at the center of a fraction of massive galaxy clusters 
(e.g. Feretti 2005 for a review). These sources are synchrotron 
emission from GeV electrons diffusing through fiG magnetic 
fields and provide the most important evidence of non thermal 
components in the ICM. 

The clusters hosting RHs are always characterized by a pecu- 
liar dynamical status indicative of very recent or ongoing merger 
events (e.g., Buote 200 1 ; Schuecker et al 200 1 ) and this suggests 
[ a connection between the mergers and the origin of non thermal 
components in the ICM. The main difficulty in understanding 
the origin of the synchrotron emitting electrons in RHs is that 
the life-time of these electrons is much shorter than the diffusion 
time necessary to cover the scale of RHs. Two main possibilities 
have been explored to explain RHs: /) the secondary electron 
models, whereby the relativistic electrons are secondary prod- 
ucts of the hadronic interactions of cosmic rays (CR) with the 
ICM (e.g., Dennison 1980; Blasi & Colafrancesco 1999), and //) 
the so-called re-acceleration models, whereby relativistic elec- 
trons injected in the intra cluster medium (ICM) are re-energized 
in situ by various mechanisms associated with the turbulence 
generated by massive merger events (e.g., Brunetti et al. 2001; 
Petrosian et al. 2001). 

The two scenarios have very different expectations for the 
basic statistical properties of RHs. In particular in the secondary 
electron models RHs should be very long-living phenomena. 
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Since the CR proton component is expected to be accumulated 
in clusters, the radio emission from secondary particles should 
be dominated at any time by the pile up of CRs protons during 
the merger history of the cluster, rather than by the latest merger 
event. On the contrary, in the re-acceleration models RHs should 
be transient phenomena with a relatively short life-time (1 Gyr 
or less) because of the finite dissipation time-scale of the tur- 
bulence. The recent discovery that galaxy clusters with similar 
X-ray luminosity have a clear bi-modal distribution in radio 
properties with only a fraction of these clusters hosting a RH 
(Brunetti et al. 2007) is in line with this latter scenario. 

Particle re-acceleration also provides a natural way to ex- 
plain the morphology and spectral properties of the few well 
studied RHs and to account for the link between RHs and clus- 
ter mergers (e.g. Petrosian 2003; Brunetti 2004; Blasi 2004; 
Feretti 2005 for reviews). The physics of collisionless turbulence 
and the process of stochastic particle re-acceleration is complex 
and rather poorly understood, however recent calculations have 
shown that there may be room for sufficient Alfvenic and magne- 
tosonic acceleration of particles in the ICM during cluster merg- 
ers (Brunetti et al. 2004; Brunetti & Lazarian 2007) thus provid- 
ing some physical support to this scenario. 

Recent calculations in the framework of the re-acceleration 
scenario have modelled the connection between RHs and cluster 
mergers in a cosmological framework, and derived the expected 
fraction of clusters with RHs as a function of the mass, redshift 
and dynamical status of the clusters (Cassano & Brunetti 2005, 
CB05; Cassano, Brunetti & Setti 2006, CBS06). These calcu- 
lations provide a unique predictive power and allow deep radio 
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observations of samples of galaxy clusters to test the scenario 
and to constrain model parameters. 

Observational studies of the statistical properties of RHs 
can be very effective in constraining the origin of RHs and, 
in general, of the non thermal components in galaxy clusters. 
Giovannini et al. (1999, GTF99) derived the occurrence of RHs 
in the XBACs sample by inspection of the NVSS (NRAO VLA 
Sky Survey, Condon et al. 1998) at 1 .4 GHz; Kempner & Sarazin 
(2001) carried out a similar study on the ACO sample by means 
of the WENSS (Westerbork Northern Sky Survey, Rengelink et 
al. 1997) at 327 MHz. Both cluster samples are complete up to 
z - 0.2. These studies show that RHs are rare at the detection 
level of the radio surveys used and that their detection rate in- 
creases with increasing the X-ray luminosity of the parent clus- 
ters: 30-35% for clusters with X-ray luminosity larger than 10"*^ 
h-^ erg/s were found to host RHs (GTF99). 

Despite these observational claims it is not clear whether the 
rarity of RHs and the increase of their occurrence with cluster 
mass is real or driven by selection biases due to the brightness 
limit of the NVSS and WENSS surveys {e.g., Rudnick et al. 
2006). 

In this paper we extend these analysis up to larger redshift, 
z < 0.4. Most important, we derive an unbiased statistics of 
RHs as we discuss the relevant radio and X-ray selection ef- 
fects and limit our analysis to a cluster sub-sample which is not 
affected by these effects. The starting sample used is this pa- 
per is the combination of the NVSS-XBAC sample at z < 0.2 
and the recent sample of 50 massive (X-ray luminous) clusters 
at 0.2 < z < 0.4 (Venturi et al. 2007) observed at 610 MHz 
with the Giant Metrewave Radio Telescope (GMRT). The sam- 
ple is presented in Sect.2 and the selection effects are discussed 
in Sect. 3. In Sect. 4 we give the results on the occurrence of RHs 
and in Sect. 5 they are compared with expectations from the re- 
acceleration scenario. Finally, in Sect.6 we discuss the main im- 
plications of this model for future low frequency radio surveys. 

A ACDM {Ho = 70 Km i ' Mpc-\ = 0.3, = 0.7) 
cosmology is adopted. 

2. Selection of the cluster samples 

In this Section we describe the sample of galaxy clusters in the 
redshift bin - 0.4 from which we start our statistical analysis 
for the presence of RH. It consists of two X-ray selected clus- 
ter samples, at low and at high redshifts, with already available 
radio follow ups for the search of RH. Being X-ray selected the 
sample may contain cooling core (CC) clusters. Because all RH 
are found in merging clusters it is well known that a possible 
anti-correlation exists between the presence of CCs and of gi- 
ant RHs at the cluster center (Edge et al. 1992; Feretti 2000). 
However, since our main aim is to perform an unbiased analysis 
of the occurrence of RH in galaxy clusters independently of their 
dynamical status, we will consider all clusters regardless of the 
presence of CCs. 

2.1. The sample at z< 0.2 

Following GTF99 at low redshifts, z < 0.2, we use the X- 
ray-brightest Abell-type clusters sample (XBACs, Ebeling et al. 
1996). The XBACs is a complete, all sky X-ray flux limited sam- 
ple of 242 clusters which are optically selected from the catalogs 
of Abell (1958) and Abell, Corwin & Olowin (1989, hearafter 
ACO) and compiled from the ROSAT All-Sky Survey. The sam- 
ple is limited to high galactic latitude (| |> 20°) and is statisti- 
cally complete for X-ray fluxes larger than 5.0x1 0" '^ergs" ' cm"^ 



in the 0.1-2.4 keV band up to a redshift 0.2, the nominal com- 
pleteness limit of the ACO clusters (Ebeling et al. 1996), but 
there are also 24 clusters at redshift larger than 0.2 which meet 
the flux criterion. 

This sample was cross-checked with the NVSS to search for 
RHs (GTF99). The NVSS is a radio survey performed at 1.4 
GHz with the Very Large Array (VLA) in configurations D and 
DnC. It has an angular resolution of 45" (HPBW), a surface 
brightness r.m.s. =0.45 mJy/beam (1 o") and covers all the sky 
north of 6 = -40°. Given the short baseline sampling available 
at the VLA, the NVSS is insensitive to diffuse sources on scales 
> 15', and this precludes the possibility to detect RH at z < 0.044 
(assuming 1 Mpc size, GTF99). Thus as a starting low redshift 
sub-sample we extract from the XBACs sample all clusters with 
0.044 < z < 0.2 and with 6 > -40. The final sample consists of 
182 clusters (excluding A 1773 and A 388 which fall in the few 
remaining gaps of the NVSS). 

2.2. The sample in the redshift range: 0.2 < z < 0.4 

The sample at larger redshift, 0.2 < z < 0.4, consists of the com- 
bination of two X-ray sub-samples with a relatively deep radio 
follow up to search for RH recently performed at 610 MHz with 
the GMRT. The two X-ray samples are the ROSAT-ESO Flux 
Limited X-ray (REFLEX) galaxy cluster catalog (Boringher et 
al. 2004) and the extended ROSAT Brightest Cluster Sample 
(eBCS) catalog (Ebeling et al. 1998, 2000). These two cata- 
logs have almost the same flux limit in the 0.1 - 2.4 keV band 
( > 3 ■ 10^'^erg s"' cm"^) and their combination yields a homo- 
geneous flux limited sample of clusters. GMRT radio follow up 
of these catalogues has been performed by our group only for 
clusters with Lx(0. 1-2.4 keV) > 5 ■ 10^"* erg s"' extracted from 
the catalogues mentioned above. Description of the two resulting 
radio-X-ray samples is given below. 

2.2.1. The REFLEX sub-sample 

The REFLEX survey covers the southern sky up to a decli- 
nation 6 - +2.5°, avoiding the Milky Way and the regions 
of the Magellanic clouds, for a total area of 13924deg2 (4.24 
sr). The sample is complete for X-ray fluxes larger than ~ 
3 ■ 10^'^erg s"' cm"^ up to z ~ 0.3; above this redshift only very 
luminous objects (with X-ray luminosities of several lO'*^ erg/s) 
have been observed (Boringher et al. 2001). In order to have a 
good u-v coverage for the radio observations at the GMRT in 
this sample we selected only clusters with 6 > -30° (Venturi et 
al. 2007). 

In Fig. [T] we report the distribution of the REFLEX clusters 
in the plane Lx - z and with red circles we highlight the clusters 
matching our luminosity and declination criteria. We obtain a 
total sample of 27 clusters. The source list is reported in Tab.[T] 
where we give (1) the REFLEX name, (2) alternative name from 
other catalogs, (3) and (4) J2000 coordinates, (5) redshift, (6) the 
X-ray luminosity in the 0. 1-2.4 keV band in unit of lO'*^ erg i"' , 
(7) information on the diffuse emission. 

Among these 27 clusters, three are clusters known to host 
RHs, i.e., A 2744, A 1300 and A 2163. From the remaining 24 
clusters in Tab. [T]we selected all clusters with no radio informa- 
tion available in the literature and we also excluded those be- 
longing to the GMRT Cluster Key Project (P.I. Kulkarni), and 
the remaining 1 8 clusters (marked with the symbol ^J in Tab. [TJ 
were observed with the GMRT (in several observational run 
from January 2005 to August 2006). The total REFLEX sample 
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REFLEX Name Alt. name RAy2ooo DECj2ooo z Lx diffuse emission 



VRXC J 0003. 1-0605 


A 2697 


00 03 11.8 


-06 05 10 


0.2320 


6.876 


none 


* RXCJ 0014.3-3023 


A 2744 


00 14 18.8 


-30 23 00 


0.3066 


12.916 


RH 


^RXCJ 0043.4-2037 


A2813 


00 43 24.4 


-20 37 17 


0.2924 


7.615 


none 


^ RXCJ 0105.5-2439 


A 141 


01 05 34.8 


-24 39 17 


0.2300 


5.762 


none 


^RXCJOllS. 1-2658 


A 2895 


01 18 11.1 


-26 58 23 


0.2275 


5.559 


none 


^RXCJOOl. 8-1336 


A 209 


01 31 53.0 


-13 36 34 


0.2060 


6.289 


RH 


^ RXCJ 0307.0-2840 


A 3088 


03 07 04.1 


-28 40 14 


0.2537 


6.953 


none 


° RXCJ 0437. 1+0043 


_ 


04 37 10.1 


+00 43 38 


0.2842 


8.989 


none 


V RXCJ 0454.1-1014 


A 521 


04 54 09.1 


-10 14 19 


0.2475 


8.178 


relic 


RXCJ 0510.7-0801 




05 10 44.7 


-08 01 06 


0.2195 


8.551 


? 


^ RXCJ 1023.8-2715 


A 3444 


10 23 50.8 


-27 15 31 


0.2542 


13.760 


core halo 


^ RXCJ 11 15.8+0129 




11 15 54.0 


+01 29 44 


0.3499 


13.579 


none 


* RXCJ 1131.9-1955 


A 1300 


11 31 56.3 


-19 55 37 


0.3075 


13.968 


RH 


RXCJ 1212.3-1816 


— 


12 12 18.9 


-18 1643 


0.2690 


6.197 


7 


^ RXCJ 1314.4-2515 


_ 


13 14 28.0 


-25 15 41 


0.2439 


10.943 


2 relics, 1 small RH 


^ RXCJ 1459.4-1811 


S780 


14 59 29.3 


-18 11 13 


0.2357 


15.531 


none 


RXCJ 1504.1-0248 


- 


15 04 07.7 


-02 48 18 


0.2153 


28.073 


7 


^ RXCJ 1512.2-2254 




15 12 12.6 


-22 54 59 


0.3152 


10.186 


none 


RXCJ 1514.9-1523 




15 14 58.0 


-15 23 10 


0.2226 


7.160 


7 


* RXCJ 1615.7-0608 


A2163 


16 15 46.9 


-06 08 45 


0.2030 


23.170 


RH 


^ RXCJ 2003.5-2323 




20 03 30.4 


-23 23 05 


0.3171 


9.248 


RH 


RXCJ 221 1.7-0350 




22 11 43.4 


-03 50 07 


0.2700 


7.418 


7 


^RXCJ 2248.5-1606 


A 2485 


22 48 32.9 


-16 06 23 


0.2472 


5.100 


none 


V RXCJ 2308. 3-0211 


A 2537 


23 08 23.2 


-02 11 31 


0.2966 


10.174 


none 


^ RXCJ 2337.6+0016 


A 2631 


23 37 40.6 


+00 16 36 


0.2779 


7.571 


none 


^ RXCJ 234 1.2-0901 


A 2645 


23 41 16.8 


-09 01 39 


0.2510 


5.789 


none 


^ RXCJ 235 1.6- 2605 


A 2667 


23 51 40.7 


-26 05 01 


0.2264 


13.651 


none 



Symbols are as follows: ^ marks the clusters observed by us with the GMRT as part of our radio halo survey (venturi et al. 2007, Venturi et 
al. in prep); * marks the clusters with extended radio emission known from the literature (A2744 Govoni et al. 2001; A 1300 Reid et al. 1999; 
A 2163 Herbig & Birkinshaw 1994 and Feretti et al. 2001). o marks the clusters without extended radio emission known from the literature 
(RXCJ0437. 1+0043 Feretti et al. 2005). The unmarked 5 clusters are part of the GMRT cluster Key Project (RI. Kulkami). 

with radio follow up thus consists of 20 galaxy clusters (here we 
exclude from this sample A 2 163 and A 209 which are already 
included in the XBACs sample, and the clusters: RXCJ 0510.7 - 
0801, RXCJ 1212.3- 1816, RXCJ 1504.1 - 0248, RXCJ 1514.9 
- 1523 and RXCJ 221 1.7 - 0350 which belong to the GMRT 
Cluster Key Project and still have no published dat^. 

2.2.2. The extended BCS sub-sample 

The 7?6>5Ar Brightest Cluster Sample (BCS; Ebeling et al. 1998) 
is a 90 per cent flux-complete sample of the 201 clusters of 
galaxies in the northen hemisphere selected from the ROSAT 
All-Sky Survey (RASS). All these clusters have flux higher than 
4.4 X 10"'^ erg cm"^ s"' in the 0.1-2.4 keV band. This sample is 
combined with a low-flux extension of the BCS (Ebeling et al. 
2000) which consists of 99 clusters of galaxies with flux higher 
than 2.8 X 10"'^ erg cm"^ s"^ in the 0.1-2.4keV band. The combi- 
nation of these two samples forms the homogeneously selected 
extended BCS (eBCS) which is statistically complete up to a 
redshiftz ~ 0.3 (Ebeling et al. 1998, 2000). 

From the eBCS catalog we selected all clusters with 15° < 
6 < 60°, and obtained a total sample of 23 clusters (Venturi et al, 
in prep). In Fig.|2]we report the distribution of the eBCS clusters 

' The exclusion of these clusters from our sample does not introduce 
any bias because they are randomly distributed in redshift and X-ray 
luminosity among the clusters in the sample 
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Fig. 1. X-ray luminosity (0.1-2.4 keV) versus z for the REFLEX 
clusters (black filled circles). Open red circles indicate the clus- 
ters belonging to our sample. 
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Table 2. Cluster sample from the eBCS catalog at z < 0.4. 



Name 




DECy2000 


z 




diffuse emission 




no 97 4.Q 8 






19 9Q 




Va^i 1 


no Of) SR 1 
uo uu Jo. i 


1 Qf, OA A\ 


u.zoou 


8 8^"^ 

O.OJ J 


none 


VAfi07 


(18 d.7 ^ 


4.Qft 90 1 9 

T^JU ^Kj iZ 


n 9890 


lU.J / 


i\ 1 1 


V22089 


09 00 45 9 


+90 55 1 3 


0.2347 


6.79 




*A773 


09 17 59.4 


+51 42 23 


0.2170 


8.097 


RH 


^A781 


09 20 23 2 


+30 26 1 5 


0.2984 


11.29 


relic ? 




no S9 SS 


1 SI S9 S9 

J 1 JZ JZ 


n 9Mn 


f\ so 

u. jy 


none 


*ZjZ001 


no AO ^7 n 


+ 1 / Uo JO 


fl Q09'> 
U. JoZJ 


1 7 70 
1 1 , ly 


complex^ 




1 n 1 7 no A 


-rjy Ul UU 


fl 9nAn 
u.zuou 


0. J7 


none 


V A 1 zL9^ 


11 ^^7 79 ^ 


+ J J 1 


n 9 1 


A 1 Q 


none 




1 J Uo UU.4 


+Z0 J\) JO 


U. jUOj 


o.yo 


none 


^A1682 


13 06 49.7 


+46 32 59 


0.2260 


7.017 


complex ' 


^Z5768 


13 11 31.5 


+22 00 05 


0.2660 


7.465 


none 


*A1758a 


13 32 32.1 


+50 30 37 


r\ f\f\ 

0.2800 


12.26 


RH 


A 1763 


13 35 17.2 


+40 59 58 


0.2279 


9.32 


none 


^Z7160 


14 57 15.2 


+22 20 30 


0.2578 


8.411 


mini halo 


^27215 


15 01 23.2 


+42 21 06 


0.2897 


7.34 


none 


^RXJ1532.9+3021 


15 32 54.2 


+30 21 11 


0.3450 


16.485 


none 


A2111 


15 39 38.3 


+34 24 21 


0.2290 


6.83 


none 


*A2219 


164021.1 


+46 41 16 


0.2281 


12.73 


RH 


A2261 


17 22 28.3 


+32 09 13 


0.2240 


11.31 


uncertain ^ 


°A2390 


21 53 34.6 


+ 17 40 11 


0.2329 


13.43 


mini halo 



Symbols are as follows: ^ marks the clusters observed with the GMRT as part of our radio halo survey (Venturi et al. 2007, Venturi et al. in prep.); 
* marks the clusters with a RH known from the literature (A 773 Govoni et al. 2001; A 1758 Giovannini et al. 2006; A2219 Bacchi et al. 2003); ° 
marks the clusters with "mini-halo" known from the literature (A 2390 Bacchi et al. 2003). For clusters A 21 1 1, A 1763 and A 2261 we analyzed 
1.4 GHz VLA archive data (Venturi et al. in prep). ' The presence of extended radio galaxies, coupled with positive residuals at the cluster centre, 
does not allow to establish the presence of possible cluster diffuse emission (Venturi et al. in prep.). ^ The presence of a RH at the cluster center is 
uncertain: only B and D-array VLA data are available and extended sources at the cluster center do not allow a firm conclusion (see discussion in 
Venturi et al. in prep.). 



in the plane Lx - z and with red circles we highlight the clusters 
which meet our selection criterifl 

Tab. 12] gives details of the clusters extracted from the eBCS 
catalogue: (1) the cluster name, (2) and (3) J2000 coordinates, 
(4) redshift, (5) the X-ray luminosity in the 0.1-2.4 keV band 
in unit of 10"*^ erg/s, (6) information on the diffuse emission. 
Among these 22 clusters (z < 0.4) 4 have already known dif- 
fuse radio emission (A 773, A 1758, A 2219, A 2390), 15 were 
observed by our group at the GMRT (September 2005-August 
2006), and for the remaining 3 clusters (A 2111, A 1763 and 
A 2261) we analyzed 1.4 GHz data from archive deep pointed 
observations at the VLA (Venturi et al. in prep.). 

The final X-ray selected sample with a radio follow up at 
z < 0.4 is thus made of 21 clusters (we also exclude A 963 which 
is already included in the XBACs sample). 

3. Radio Selection effects 

The final sample of X-ray selected clusters with a radio follow 
up from which we start our analysis is made by the combina- 
tion of the XBAC-NVSS, the REFLEX-GMRT and the eBCS- 
GMRT sub-samples and consists of 220 galaxy clusters in the 
redshift interval 0.044 < z < 0.4. The distribution of the clusters 
in the - z plane is reported in Figl3] where we also highlight 
with red open circles the clusters with giant RH (GTF99, Venturi 
et al. 2007, Venturi et al. in prep.). 

^ Formally the eBCS sub-sample contains also RXCJ2228. 6+2037 
which is at z > 0.4 and is not included in the following analysis. 



In order to perform a statistical census of the fraction of clus- 
ters with giant RHs in different ranges of X-ray luminosity and 
redshift it is crucial to study the selection effects on the detec- 
tion of RHs. These biases come from the strategy adopted in the 
radio follow up of the different X-ray sub-samples. 

In the following we discuss the selection effects on the NVSS 
and GMRT cluster samples and derive a radio - X-ray sub- 
sample of clusters suitable for an unbiased statistical study of 
RHs. 



3. 1 . The XBAC-NVSS sub-sample 

As previously reported, the NVSS is surface brightness-limited 
and RHs with a lower brightness are lost. 

There is a clear evidence of the existence of a correlation 
between the synchrotron radio power at 1.4 GHz (f 1.4) of ha- 
los and the X-ray luminosity (Lx) of the parent clusters (Liang 
et al. 2000; Feretti 2000, 2003; EnBUn and Rottgering 2002; 
Cassano et al. 2006). This correlation is not affected by selec- 
tion effects, at least for high X-ray luminosities (Log(Lx) > 44.6 
erg/s;e.g., Clarke 2005, Dolag 2006, Brunetti et al. 2007) and 
it implies that the average brightness of RHs increases with the 
X-ray luminosity of the parent cluster (Feretti 2005; CBS06). If 
we scale the radio brightness with the cluster X-ray luminosity, 
the sensitivity limit in radio surveys will set the minimum sur- 
face brightness of a RH to be detected, and then a limiting X-ray 
luminosity of the hosting clusters. 
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Fig. 3. X-ray luminosity (0.1 - 2.4 keV) versus z for the total sample of clusters (XBACs-i-REFLEXH-eBCS). Open red circles 
indicate the clusters with known giant RHs. The lines give the lower limit to the cluster X-ray luminosities in the case a) (black 
line), in the case b) (dotted line), and in the case c) (blue dot-dashed Une; 2cr, blue dashed line; 3cr) as described in the text. 
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Fig. 2. X-ray luminosity (0. 1-2.4 keV) versus z for the REFLEX 
clusters (black filled circles). Open red circles indicate the clus- 
ters belonging to our sample. 



An additional argument that should be used to discuss the 
selection of RHs in brightness limited surveys is their bright- 
ness distribution. In general, RHs are characterized by brightness 
profiles which smoothly decrease with distance from the cluster 



center (e.g., Govoni et al. 2001). The outermost low brightness 
regions of RHs are lost in shallow surveys, however what is im- 
portant here is the capability of these surveys to detect the cen- 
tral, brightest part of the diffuse emission. Indeed the detection 
of these regions allows to claim the presence of diffuse emis- 
sion from clusters and to select samples of RH candidates for 
a deeper radio follow up. By making use of several well stud- 
ied RHs, it was shown (Brunetti et al. 2007) that their integrated 
brightness profiles are quite similar, provided that their radial 
distance is normalised to the size of the different RHs. In par- 
ticular it has been shown that about half of the total radio flux 
density measured in RHs is contained in about half radius and 
this information can be used in discussing the selection of RHs 
in brightness limited surveys. 

Thus in discussing selection biases from the NVSS we use 
both the Pi 4-Lx correlation and the constraints from the RH 
brightness profiles, and consider three possible approaches: 

- case a): we consider the case of RHs with a fixed radius, 
Rh - 500 kpc. By making use of the Pi ^-Lx correlation, for 
each redshift we calculate the minimum X-ray luminosity of 
a cluster that can host RHs with an expected average radio 
brightness of 0.45 mJy/beam within half radius. The result- 
ing minimum X-ray luminosity is reported as a function of 
redshift in Figl3](solid black line). 

- case by. we consider an additional correlation which has 
been recently found between the size of RHs and the syn- 
chrotron radio power at 1.4 GHz, f i 4 oc 7;^i8±o.68 (c^ssano 
et al. 2007, C07 hearafter): clusters with low X-ray luminosi- 
ties should host smaller RHs than those with higher X-ray lu- 
minosities. This makes the detection of RHs in clusters with 
lower X-ray luminosity easier with respect to the previous 
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Fig. 4. Example of injection of fake RHs. We assume a size of 500 kpc and a flux of 0, 28, 32, 36 and 45 mJy (from top left to bottom 
right panel). Ther.m.s. in the images is of 0.45 mJy/beam, the HPWB is 45"x45", contours are 1. Ox (-1, 1, 1.4,2,4,8,32, 125,500) 
mJy/beam. 



case a). For each z, we compute the minimum radio power 
of RHs with an average radio brightness within half radius 
(which depends on the radio power) of 0.45 mJy/beam. The 
corresponding minimum X-ray luminosity of the clusters 
which may host these RHs is reported in Fig|3]as a function 
of redshift (red dotted curve). 
- case c): we consider RHs with fixed size, Rh = 500 kpc, and 
with a typical brightness profile taken from Brunetti et al. 
(2007). For each redshift, we integrate the average brightness 
over an area as large as 4 times the NVSS beam and calculate 
the corresponding 1 .4 GHz power of RHs whose brightness 
equals 2 and 3 times the rms of the NVSS (0.9 and 1.35 
mJy/beam). 

The corresponding minimum X-ray luminosity (from the 
PiA-Lx correlation) of clusters which may host these RHs is 
reported in Fig. [3] as a function of redshift (blue dot-dashed 
and blue dashed lines mark the case of 2 and 3 times the 
NVSS-rms, respectively). 

In all the cases the NVSS is efficient in detecting RHs in 
clusters with Log(Lx) > 44.6 erg/s. 

In order to strengthen the reliability of the above considera- 
tions on the efficiency of the NVSS in catching RHs, we use a 
more direct approach. Adopting the same procedure of Brunetti 
et al. (2007) we choose a cluster without RH and inject a "fake" 
RH in a region close to the pointing and without any significant 
contribution from real sources (upper left panel in Fig. |4]l. The 



"fake" RH was injected in the u-v data by means of the task 
UVMOD in AIPS and the resulting new data set was imaged us- 
ing standard procedures. The brightness profile of the fake RH 
was modelled with a set of optically thin spheres with varying 
radius and flux density (Brunetti et al. 2007; Venturi et al. in 
prep). The u-v data-set used for this procedure was chosen from 
VLA archive observations carried out in the D configuration and 
it was edited so as to mimic a typical NVSS field in terms of u-v 
coverage, sensitivity and observing beam. 



In Fig|4]we report the results from the injection of fake RHs 
with fixed radius, Rh = 500 kpc and different total flux densi- 
ties, and assuming z=0.1 (which is a mean value for the NVSS 
sample). From left top panel to right bottom panel we report the 
results from injections of fake RHs with increasing total flux 
density. We find that diffuse emission is unambiguously imaged 
for RHs with flux density between 28 and 32 mJy. At the as- 
sumed redshift this corresponds to 1 .4 GHz luminosity of RHs 
of f 1 4 ~ 7-8-10^^ Watt/Hz, which corresponds to cluster X-ray 
luminosity Lx ~ 4- ■ 10^"^ erg/s (once the f 1 4 - Lx is assumed), 
in good agreement with our estimates of the RH-selection ef- 
fects in the NVSS in the most conservative case a) (Rh - 500 
kpc, black line in Fig.O, which therefore will be adopted in the 
following of the paper. 
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3.2. The GMRT sub-sample 



The case of the GMRT sub-sample is different from that of the 
NVSS sub-sample, because the observations carried out by our 
group (Venturi et al. 2007, Venturi et al in prep) are much deeper 
(Icr ~ 35 - 100 /iJy beam"') and aimed at the detection of RHs. 
This guarantees that the detection of extended diffuse emission 
at the level expected from the f i 4 - Lx correlation is not biased 
by sensitivity limit of those observations. 

Brunetti et al. (2007) derived solid upper limits to the radio 
powers of clusters (in the GMRT sample) without evidence of 
diffuse radio emission. These upper limits lie about one order 
of magnitude below the RH luminosities expected on the basis 
of the Pi 4 - Lx correlation. This allows to conclude that the 
radio follow up performed for the GMRT sub-sample does not 
introduce relevant biases to the detection of RHs in the cluster 
sample. 



4. Results: fraction of galaxy clusters with radio 
halos 

Based on the results reported in the previous Sections, we define 
a new cluster sample which is made by the GMRT sub-sample 
(Sect. 2.2) and by all clusters of the XBACS-NVSS sub-sample 
(Sec. 2.1) not affected by biases in the detection of RHs, namely 
all those above the curves in Fig. [3] 

In Fig.[3]the clusters hosting a giant RH are represented with 
open red circles: it is clear that RHs most likely occur in high 
X-ray luminosity clusters, and this is consistent with previous 
claims (GTF99). 

To obtain a statistically meaningful estimate of this "be- 
haviour" we focus on the most conservative sub-sample, which 
is made by all clusters which lie above the upper curve in Fig.[3j 
this curve is derived by adopting the case a) (Sec. 3.1), that is 
supported by UVMOD simulations (Fig.|4|. 

As already anticipated the REFLEX and eBCS samples are 
complete at the X-ray flux limit up to z a: 0.3, and Fig. [3] 
shows that incompleteness strongly affects the population of 
clusters with X-ray luminosity typical of the GMRT sub-sample 
at z > 0.32. Thus we limit our analysis to the redshift range 
0.044 - 0.32 and calculate the fraction of clusters with RHs in 
the population of low luminosity (LL) and high luminosity (HL) 
clusters. We chose Lx = lO"*"*'^ erg/s as transition value between 
the LL and HL samples. Such value provide good statistics in 
both samples and ensure that the HL clusters cover the lumi- 
nosity interval where the occuiTence of RHs seems to increase 
(Fig. [3]). The LL sub-sample is thus composed by those clusters 
with X-ray luminosity above the black curve in Fig.[3]and below 
10'*'* '^ erg/s. 

We found that the fractions of clusters with RHs are fun - 
0.406 + 0.1 12 and /rh ^ 0.075 + 0.038 in the HL and LL sub- 
samples, respectively (Icr Poissonian eiTors). 

This is the most conservative case, indeed it is worth men- 
tioning that by considering LL sub-samples defined by other 
selection-curves (case b) and c)), smaller values of /rh are ob- 
tained in the LL clusters: ^ 0.063 ± 0.026 in the case b) (red 
dotted curve in Fig. [3]l and fun - 0.074 ± 0.030 in the case c) 
(Isigma, blue dot-dashed line in Fig.|3]l. 

In all these cases there is a clear evidence of an increase of 
the fraction of clusters with RHs with increasing the X-ray lu- 
minosity. To test the strength of this result we consider the most 
conservative case and run Monte Carlo calculations. Specifically 
in this case the sample is composed by 85 clusters (32 are HL 
clusters and 53 are LL clusters) with 17 RHs, of which 13 are 



KVSS+GKRT 
10= trails 




Nji in the upper bm 



Fig. 5. Montecarlo realization of 10^ trials in the HL clusters in 
flie NVSS-hGMRT sample. 



hosted by HL clusters and 4 by LL clusters. We randomly assign 
17 RHs among the 85 clusters of our sample and count the num- 
ber of RHs which fall into the HL sub-sample in our Monte Carlo 
trials. In Fig. |5] we report the distribution of the number of RHs 
in the HL sub-sample obtained after 10^ Monte Carlo trials. The 
distribution can be fitted with a Gaussian function with a central 
value of ~ 6.29 and a standard deviation of ~ 1.81. This means 
that the observed value of 13 RHs in the HL clusters is at 3.7cr 
from the value expected in the case of RH occurrence indepen- 
dent of cluster X-ray luminosity, and thus that the probability 
to have the observed fraction of clusters with RHs by chance is 
< 0.2 ■ 10"^^. Here we have assumed the presence of a RH in the 
cluster A 2261. This does not affect our result as the observed 
jump remains highly significant (~ 3.4cr) even by assuming that 
A 2261 does not have a RH. 

This shows that our result is statistically significant and sug- 
gests, for the first time, that there is a real (i.e., physical) increase 
of the occurrence of RHs with increasing the X-ray luminosity 
of the hosting clusters in an X-ray selected sample. 

In principle, because we are using X-ray selected cluster 
samples, it is not immediate to conclude if the increase of the 
fraction of clusters with RHs across the LL and HL clusters is 
purely driven by the X-ray luminosity of clusters or is due to 
possible difference in the redshift distribution of clusters in the 
two HL and LL sub-samples. In Fig|6]we report the redshift dis- 
tributions of clusters in the two sub-samples. HL clusters have 
typically larger redshift because they can be detected by X-ray 
surveys at larger cosmic distance. More specifically, the mean 
redshift of LL clusters is ~ 0.16, while that of HL clusters is 
~ 0.23. However the difference in redshift distribution is rel- 
atively small (~ 0.74 Gyr in cosmic time) and thus it is very 
unlikely that the increase (a factor of ~ 5) of the fraction of clus- 
ters with RHs observed in the two sub-sample is due to some 
cosmic evolution of the population of RH sources. In addition, 
regardless of the origin of the particles emitting RHs, it should 
be mentioned that the radiative life-time of these particles de- 
creases with redshift due to IC losses and thus, in general, the 
fraction of clusters with RHs would be expected to decrease 
with cosmic look back time. Finally we note that the increase 
of the fraction of clusters with RH with cluster X-ray luminosity 
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appears in Fig. [3] also considering the GMRT sample alone, in 
which case clusters are essentially at the same redshift. Thus, al- 
though it is true that both luminosity and redshift may influence 
the occurrence of RHs (and model calculations that account for 
both effects are discussed in Sec. 5), we believe that the X-ray 
luminosity plays the major role in our cluster sample. 

One possibility to further test this issue and to minimize the 
possible effect of the redshift on the observed fraction of clus- 
ters with RHs is to consider the clusters of the XBAC-NVSS 
only which lie above the most conservative radio selection-curve 
in Fig. [3j this sample is selected within 0.044 < z < 0.2. In 
this case, however, it should be considered that the smaller clus- 
ter statistics and range of X-ray luminosities in the sample are 
expected to decrease the statistical significance of any possible 
trend observed for the fraction of RH clusters with cluster prop- 
erties. Also in this case we derive the fraction of clusters with 
RHs in LL and HL clusters but to maintain a sufficient cluster 
statistics in both the sub-samples we adopt an X-ray luminosity 
threshold of - 10'*'* l 

We find /r/^ ^ 0.43 + 0.17 and /rh ^ 0.088 + 0.051 in the HL 
and LL clusters, respectively, which again suggests an increasing 
fraction of clusters with RHs with increasing X-ray luminosity. 

As in the previous case, in order to test the statistical strength 
of this results, we perform Monte Carlo calculations by assum- 
ing equal probability to have RHs in the clusters of the sample. 
In this case we have a sample of 48 clusters (14 are HL clusters 
and 34 are LL clusters) with 9 RHs, 6 of these RHs are in the 
HL sub-sample. In Fig|7]we report the distribution of RHs in the 
HL sub-sample as counted after 10^ trials. The distribution is fit- 
ted with a Gaussian function which peaks at ~ 2.49 and with a 
standard deviation of ~ 1 .29, thus the observed value of 6 RHs 
is at ~ 2.7cr from that expected in the case of equal probability 
to have RHs with cluster X-ray luminosity; this means that the 
probability to obtain by chance the observed fraction of clusters 
with RHs is < 9 ■ 10 I 

Thus, although the significance of the result is reduced with 
respect to the case in which we consider the full NVSS+GMRT 
sample, the observed jump in the occurrence of RHs with the 
X-ray luminosity is still statistically meaningful and further sug- 



gests that the effect is physically driven by cluster X-ray lumi- 
nosity (or mass). 

4. 1 . The effect of cooling core clusters on the observed 
statistics 

It should be mentioned that the Pi 4 - Lx correlation we used 
in Sec. 3 to derive the capability of the NVSS in catching RHs, 
is found for RH clusters which do not have a CC. In principle, 
the X-ray luminosity of these clusters should be corrected for 
the contribution of the CC when using the f j 4 - Lx correlation. 
Due to the lack of an adequate X-ray follow up of the clusters in 
our sample, it is impossible to identify those affected by a CC. 
However the fraction of CC clusters is found at a level of ~ 30% 
for BCS clusters with > 4 ■ lO'*'* at z ~ 0. 15 - 0.4 (Bauer et al. 
2005) and this is unlikely to affect our main results. Indeed once 
corrected from the CC contributions some HL clusters may be- 
come LL clusters and some LL clusters may fall below the X-ray 
luminosity threshold and leave the selected sample. Because the 
LL and HL clusters are comparable in number, the occurrence 
of RH in both HL and LL clusters is expected to remain un- 
changed (although the fraction of RH in HL clusters may slightly 
increase). In order to better quantify this effect we work out 
Montecarlo simulations. We assume that 30% of galaxy clusters 
of our sample are CCs and randomly assign CCs among clus- 
ters without RHs. Once the CC is assigned to a cluster we apply 
a correction to its X-ray luminosity which is reduced by a fac- 
tor randomly chosen between 1 .4 and 2.5 (consistent with recent 
observations, e.g., Zhang et al. 2006, 2007; Chen et al. 2007) and 
re-evaluate the statistics of RH in the "corrected" LL and HL 
samples. In Fig. [8] a) we report the distributions of the fractions 
of clusters with RH in the "corrected" sample, LL (left distribu- 
tion) and HL (right distribution), obtained after having applied 
our procedure in 10^ Montecarlo trials. The fractions of clusters 
with RHs in both LL and HL clusters only slightly increase with 
respect to the previous statistical calculation (Sect.4). In Fig. |8] 
b) we report the distribution of the ratio between the fraction of 
clusters with a RH in the "corrected" HL and LL samples ob- 
tained from 10^ Montecarlo trials: it is clear that the bulk of the 
values is consistent with that from the previous statistical anal- 
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Fig. 8. a) Distributions of the fractions of clusters with RH in LL (left histogram) and HL (right histogram) sub-sample, obtained 
after 10^ Montecarlo trials. The red filled circles rapresent the corresponding value of /rh obtained without CC correction, b) 
Distribution of the ratio between the fraction of clusters with RH in HL and LL clusters obtained from 10* Montecarlo trials, the 
red filled circle being the value obtained without CC correction. 
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Fig. 9. Expected fraction of clusters with RHs in the HL and LL 
sub-samples (different colored lines) vs the observed occurrence 
of RH with Icr errors (shadowed regions). 

ysis. These results, based on viable assumptions, show that the 
main finding of our work, the increase of the RH occurrence with 
the X-ray luminosity of galaxy clusters, might not be affected by 
the possible presence of CC clusters in the sample. 

5. A comparison between model expectations and 
observations 

In this Section we provide a basic comparison between present 
observations and the expectations of the re-acceleration model 
in its simplest form. In the previous Section we showed that the 



fraction of clusters with RHs increases with the cluster X-ray 
luminosity and this provides an important piece of information 
that should be explained by the theoretical pictures for the ori- 
gin of RHs. CB05 modelled the statistical properties of RHs in 
the framework of the merger-induced in situ particle accelera- 
tion scenario. By adopting the semi-analytic Press & Schechter 
(1974) theory to follow the cosmic evolution and formation of a 
large synthetic population of galaxy clusters, it was assumed that 
the energy injected in the form of magnetosonic waves during 
merging events in clusters is a fraction, 77,, of the PdV work done 
by the infalling subclusters in passing through the most massive 
one. Then the processes of stochastic acceleration of the rela- 
tivistic electrons by these waves, and the ensuing synchrotron 
emission properties, have been calculated under the assumption 
of a constant magnetic field intensity averaged within a 1 Mpc^ 
volume. CBS06 have extended this analysis, by including a scal- 
ing of the magnetic field strength with cluster mass in the form 
B oc M*. They showed that the observed correlations between 
the synchrotron radio power of a sample of 17 RHs and the X- 
ray properties of the hosting clusters can be reproduced in the 
framework of the re-acceleration model for b > 0.5 and typical 
juG strengths of the average B intensity. Those values provide 
a working framework in the range of (B, b) model parameters 
under which the statistical properties of RHs, namely the occur- 
rence of RHs with cluster mass and z, the luminosity functions 
and number counts, have been extensively calculated. 

An important finding of this work was that generally the 
probability to form giant RHs increases with cluster mass, which 
is a unique expectation of the re-acceleration scenario. In partic- 
ular an abrupt increase of the fraction of clusters with RHs is ex- 
pected across cluster masses of ~ 2 ■ lO''' M©, since the turbulent 
energy injected during cluster mergers scales with the cluster 
thermal energy (which roughly scales as ~ M^^^), and turbu- 
lence in more massive clusters is injected on larger volumes (see 
CB05). 

Specifically, the most important goal of this Section is to test 
whether this re-acceleration scenario can potentially match the 
increase of the fraction of clusters with RHs as derived from 
our cluster sample. We therefore adopt the same approach out- 
lined in CB05 and CBS06 and compare model expectations with 



10 



R. Cassano et al.: Revised statistics of radio halos and the re-acceleration model 



the results obtained in the previous section (case NVSS+GMRT 
samples). 

In order to have a prompt comparison between model expec- 
tations and the data we calculate the expected fraction of clusters 
with RHs in different redshift bins by assuming the same redshift 
distribution of galaxy clusters as in Fig. |6] It should be men- 
tioned that the model use virial masses (M,,) and thus to compare 
our expectations with observations we use the observed Lx - M,, 
COIT elation (Eq. 7 in CBS06): Lx oc mI-*\ This coiTelation is 
used to convert the range of X-ray luminosity of LL and HL 
clusters into mass ranges in which to compute the expected frac- 
tion of clusters with RHs. 

In Fig. |9] we report the derived for different values of 
the model parameters (different lines) as a function of the X-ray 
luminosity overlaid on the observed values of /rh (shadowed 
region). The model parameters are the magnetic field intensity, 
B<M>, of a cluster with virial mass < M >= 1 .6 ■ 10'^ M© and the 
slope b of the scaling of B with M,, (see CBS06 for details). In all 
the three cases a value of 77, ^ 0.18 is adopted: this value falls in 
the range of 77, that has been shown to allow the re-acceleration 
model, in the form developed by CB05 and CBS06, to be con- 
sistent with the production of RHs in ^ 1/3 of massive galaxy 
clusters. We find that the expected behaviour of the occuiTence 
of RHs with cluster mass is consistent with the observed one for 
all the considered values of parameters. 

It should be mentioned that the use of the Lx - My corre- 
lation implies that Lx reported in Fig|9]is a good proxy for the 
mass. This Lx-My correlation is obtained in CBS06 by using the 
HIFLUGCS cluster sample (Reiprich & Bohringer 2002) which 
contains CC and non CC clusters and which is characterized by 
a relatively large scatter in the mass, i.e. ^45% (Reiprich & 
Bohringer 2002). Very recent results on the HIFLUGCS sample 
showed that most of this scatter is due to the mixing of CC and 
non CC clusters in the sample, with CC and non CC clusters 
following correlations with compatible slopes (within 1 cr) but 
different normalization^ (Chen et al. 2007; Reiprich & Hudson 
2006). 

Thus the use of the Lx - My coiTelation of CBS06 implies that 
the masses of non CC clusters and CC clusters with a fixed Lx 
would be statistically under-estimated and over-estimated, re- 
spectively, and this would provide a compromise in the calcula- 
tion of the masses spanned by the LL and HL clusters. A more 
detailed comparison would require the knowledge of the virial 
masses for all the clusters of our sample, which is however not 
available at the moment. Nevertheless this is unlikely to affect 
our results significantly. Indeed one simple check is to use the 
/rh in Fig. [8](/.e., obtained from Montecarlo simulations after 
statistical correction of the Lx of CC clusters) and to apply a 
systematic positive correction of « 30 % to the mass of LL and 
HL clusters with respect to the case of the Lx - My correlation 
ofCBS06Q.Inthis case we checked that the jump in /rh and the 
values of /rh themselves in Fig. [8] can be still well reconciled 
with the model expectation^ 



' This is shown for both the Lx - Msoo and Lx - M200 correlations. 
This correction has been estimated from the Lx - M200 correlation 
for non CC clusters in Fig. 2 of Reiprich & Hudson 2006. 

^ We note that in this case, because the cluster masses (energetics) are 
shghtly larger, the model formally reproduces /rh for a slightly smaller 
value of T],, Tj, x 0.16. 



6. Radio Halos in future low-frequency radio 
surveys 

In the previous Section we showed that the re-acceleration sce- 
nario may explain the observed increase of the fraction of clus- 
ters with RHs with cluster X-ray luminosity (mass). On the other 
hand, from Fig|9] it is evident that with the present surveys it 
is difficult to make a more quantitative statement, i.e., the poor 
statistics does not allow to disentangle among different configu- 
rations of the model parameters. 

In a few years LOFAR and LWA will survey galaxy clus- 
ters with unprecedented sensitivity, thus increasing the obser- 
vational statistics and hopefully providing a much clear situa- 
tion. On the other hand these radio telescopes will observe at 
frequencies ^ 10 times lower than the present facilities and this 
could make the interpretation of the observations more tricky. 
Indeed another natural expectation of the re-acceleration sce- 
nario is that the fraction of clusters with RHs depends on the ob- 
serving frequency. In particular, it is expected that the number of 
RHs should increase at low radio frequencies, a; 100 MHz, since 
low energy electrons (with Lorentz factor y » 3000) are able to 
radiate synchrotron emission at these frequencies and may be 
re-accelerated also in less energetic merger events (CBS06). For 
energetic reasons it is also expected that this increase should be 
more pronounced in the case of less massive clusters which thus 
should host a sizable fraction of the population of RHs emitting 
at low radio frequencies. 

The new generation of low frequency interferometers, 
LOFAR and LWA, represents therefore a unique possibility to 
discover a new class of Radio Halos predicted on the basis of the 
re-acceleration scenario which may emerges only at low radio 
frequencies. In particular, it is expected (CBS06) that the num- 
ber of low frequency RHs should be « 10 times that of "classical" 
RHs (those emitting at siGHz frequencies), thus the discovery of 
an excess of RHs at low frequencies with respect to the extrap- 
olation of the number counts of present RHs would be a prompt 
confirmation of this scenario. 

In FigHO] we report some preliminary calculations (based 
on CBS06, Sect. 7) of the probability to have RHs at different 
radio frequencies (1.4 GHz, 240 MHz and 150 MHz, see fig- 
ure caption). We found a large increase of the expected fraction 
of clusters with RHs at low radio frequencies (150-240 MHz); 
this increase is even more striking for less massive clusters, 
i.e., M < 10'^ Mq. Indeed if the fraction of clusters with RHs 
in the redshift bin 0-0.1 increases between 1400 - 150 MHz by 
a factor of ~ 2 for M ~ 2 - 3 ■ 10'^ Mq, this increase is a fac- 
tor of > 10 for M < 10'^ Mq. Furthermore, at higher redshifts, 
z ~ 0.5 - 0.6, this increase is even larger, a factor of ~ 4 and 
> 10^ in the case of M ~ 2 - 3 ■ 10'^ Mo and M < lO^^ Mq, 
respectively. From these results, we reach the important conclu- 
sion that the increase of the fraction of clusters with RHs with in- 
creasing cluster mass (or X-ray luminosity) should become less 
striking at low radio frequencies. In particular, considering the 
mass range in which present radio surveys have observed RHs, 
i.e., IO'^Mq < M < 4 ■ IO'^Mq, and considering two population 
of clusters with a separation mass ~ 2 ■ 10'^ Mq, the jump in the 
expected fraction is a factor of ~ 4 at 1 .4 GHz and is reduced to 
a factor of ~ 2 at 150 MHz. Thus, in general, LOFAR and LWA 
surveys are expected to find a smaller increase of the fraction of 
clusters with RHs with the cluster mass with respect to present 
surveys. 
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Fig. 10. Expected fraction of clusters with RHs at three different radio frequencies: 1.4 GHz (blue lines), 240 MHz (red lines) and 
150 MHz (black lines) in the redshift bin z ~ - 0.1 (left panel) and z ~ 0.5 - 0.6 (right panel). Calculations have been performed 
assuming: b = 1.5, B<m> = 1.9 yuG and t], = 0.18. 



7. Conclusions 

The statistical properties of RHs are an important piece of infor- 
mation for the present understanding of non-thermal phenomena 
in galaxy clusters: do all clusters have a RH? Does the fraction 
of clusters with RHs depend on cluster mass and redshift? 

In the past a few seminal attempts were made to address 
these points (Giovannini et al. 1999; Kempner & Sarazin 2001), 
however it was not clear how much those results were affected 
by the sensitivity limit of the adopted radio surveys. In partic- 
ular it was unclear whether the absence of RHs in low X-ray 
luminosity clusters could be ascribed to the steepness of the ra- 
dio power-X-ray luminosity correlation of RHs combined with 
the sensitivity limit of the adopted surveys (Kempner & Sarazin 
2001). 

In this work we presented an unbiased statistical analysis 
of RHs in a large sample of clusters with radio information. 
We combined the sub-sample of XBACs clusters at z < 0.2, 
checked for the presence of diffuse radio emission in the NVSS 
by GTF99, with two sub-samples of clusters in the redshift range 
0.2 - 0.4 (exti-acted from the REFLEX and eBCS catalogs) in- 
spected at 610 MHz and at 1.4 GHz with deep GMRT and 
VLA observations, the GMRT sub-sample (Venturi et al. 2007; 
Venturi et al. in prep.). 

Brunetti et al. (2007) showed that the bulk of clusters of the 
GMRT sample does not show any hint of Mpc scale diffuse radio 
emission at their centre at the level of presently known RHs. For 
the clusters without extended radio emission upper limits to their 
radio power were obtained: they are about one order of magni- 
tude below the values expected on the basis of the radio power- 
X-ray luminosity correlation, and this allowed to conclude that 
such correlation is real, i.e. not driven by observational biases at 
least for Lx > 5 ■ 10'^'* erg/s (the lower X-ray luminosity of the 
GMRT sample). 

Assuming that RHs follow this correlation, in this paper we 
have obtained the minimum X-ray luminosity of clusters which 
may host a RH detectable in the NVSS as a function of z. We 



followed a number of approaches based on the brightness pro- 
files of well studied RHs. In the most conservative case (case 
a). Sect. 3.1) we successfully test this approach by means of in- 
jection of fake RHs in a NVSS-like UV dataset. Based on these 
results for the NVSS and on those of Brunetti et al. (2007) for 
the GMRT sample, we were able to select a sample of clusters 
(NVSS+GMRT) which is not affected by observational biases 
and which is suitable for a statistical analysis. 

The main result of our analysis of this sample is that we find 
an increase of the fraction of clusters with RHs with the cluster 
X-ray luminosity. More specifically, we find that the percentage 
of clusters at z * 0.044 - 0.32 with 4 ■ IQ'^'^ergls < Lx < 8 ■ 
10"^ ergi s hosting RHs is 7.5 + 3.8%, while that of clusters with 

> 8 ■ 10'*'* erg/s is 40.6 + 1 1 .2%. We tested the significance of 
this result by means of Montecarlo trials which allow us to con- 
clude that the observed jump is real, with a significance of 3.7cr. 
In Sec. 4.1 we also showed that this result is not appreciably 
affected by the possible presence of CC clusters in the sample. 

We showed that the increase of the fraction of clusters with 
RHs can be understood in the framework of the re-acceleration 
scenario. This increase is a natural expectation of this scenario 
since a more efficient particle acceleration is triggered in massive 
(and X-ray luminous) clusters during merger events. 

Since this scenario is at present in line with several observa- 
tions, it is extremely intriguing, in view of forthcoming low fre- 
quency interferometers, such as LOFAR and LWA, to investigate 
the model expectations at low radio frequencies, i.e. few hun- 
dreds of MHz. In Sec. 6 we showed that this model expects that 
the population of low frequency RHs is 10 times larger than 
that of "classical" RHs (those emitting at « GHz frequencies) 
and we showed that the increase of the probability to have RHs 
emitting at low radio frequencies is maximized in the case of 
clusters with lower masses (M < 10'^ Mq). This implies that the 
increase of the fraction of clusters with RHs with cluster mass 
as measured by LOFAR and LWA surveys is expected to be less 
pronounced than that observed by present surveys. 
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